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Abstract 



I In the framework of perturbative QCD and potential models of 

heavy quarkonium, total cross-sections and differential characteristics 
of hadronic production for different spin states of the Sc-meson family 
are calculated at energies of LHC. Theoretical predictions for masses 
^ and branching ratios of Be decays are given, which allows one to make 

■ estimates for an expected number of reconstructed events with Be- 

en 

o 

: Introduction 

p • Along with a possible direct observation of higgs particles and new particle 



interactions beyond the Standard Model, a detail investigation of processes 



D . with 6-quarks, involving mechanisms caused by high orders of the QCD per- 

^ ' turbation theory, is one of the most important tasks of LHC physics pro- 

gram. Certainly, in the hadron decays with 6-quarks the precision study of 
^ ■ electroweak model parameters and the observation of "new" physics effects 

■ are accessible. The latters appear as virtual corrections to the leading ap- 

proximation and, thus, must be reliably separated from the QCD high order 
corrections. So, the investigation of CP-violation in the heavy quark sector 
supposes an observation of more than 10^ decays of B mesons. While per- 
forming this experiment it is also possible to solve another problem such as 
the observation of meson and investigation of its properties, which are 
attracting a significant interest |jl|. With respect to the cross-section for the 
production of B mesons, the yield of Be is suppressed by a hard production 
of the additional heavy c-quark as well as a small probability of the (be) state 
formation, so that with a{Bf)/a{bb) ~ 10^'^ one can expect about 10^ events 
with Be per 10^ i?-mesons. 

The (fee) system is a mixed flavour heavy quarkonium, whose states have 
no annihilation decay modes due to the electromagnetic and strong interac- 
tions. Therefore, the excited levels of (6c), lying below the threshold of decay 
into the pair of heavy B and D mesons, follow the radiation decays into the 
pseudoscalar ground state, B'^ meson, which is a longlived particle, decaying 
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due to the weak interaction, t{Bc) — 0.55 ± 0.15 psec. A nonrelativistic 
motion of heavy quarks in the (be) system causes a rehable apphcation of 
phenomenological potential models worked out for the charmonium (cc) and 
bottomonium (bb) , where the precision of level mass prediction reaches 30 eV 
[|, Weak exclusive decays of are described in the framework of quark 
models for form-factors of different currents ^, which allows one also to 
determine the total width of this state in the form of a sum over all possible 
decay modes calculated. The value obtained in this way agrees with esti- 
mates of inclusive Be decays described in the framework of the quark-hadron 
duality approach and operator product expansion 0. In this consideration 
one takes into account the virtualities of heavy quarks in the initial state, 
annihilation decay channels and Pauli interference effects for the c-quark, 
being a product of the 6-quark decay mode, with the c-quark from the initial 
(be) state. 

The Be production mechanisms can be reliably investigated in the frame- 
work of perturbative QCD for a hard production of two heavy quark pairs 
and a potential model of a soft nonperturbative binding of b and c quarks 
with a small relative momentum in the system of quark mass centre. In 
e"'"e~-annihilation at large energies (M^^/s <^ 1), the consideration of lead- 
ing diagrams for the _Bc-meson production gives the factorized scaling result 
for the differential cross-section over the energy fraction carried out by me- 
son, da/dz = a(bb)-D{z), where z = 2EB^/y/s 10) @- This distribution allows 
the interpretation as the hard production of heavier 6-quark with the sub- 
sequent fragmentation into B^, so that D{z) is the fragmentation function. 
In this approach one can obtain analytic expressions for the fragmentation 
functions into the different spin states of S, P, and D-wave levels In 
photon-photon interactions for the Be production in the leading approxima- 
tion of perturbation theory 0, one can isolate three gauge-invariant groups 
of diagrams, which can be interpreted as 

1) the hard photon-photon production of bb with the subsequent frag- 
mentation of 6 — > B^{nL), where n is the principal quantum number of the 
(6c)-quarkonium, L = 0, 1 ... is the orbital angular momentum, 

2) the corresponding production and fragmentation for the c-quarks, and 

3) the recombination diagrams of (6c)-pair into B^, wherein the quarks 
of different flavours are connected to the different photon lines. 

In this case, the results of calculation for the complete set of diagrams in 
the leading order of perturbation theory show that the group of 6- fragmentation 
diagrams at high transverse momenta pt(-Bc) ^ can be described by 
the fragmentation model with the fragmentation function Dj^^^+{z), calcu- 
lated in the e"'"e~-annihilation. The set of c-fragmentation diagrams does not 
allow the description in the framework of fragmentation model. The recom- 
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bination diagrams give the dominant contribution to the total cross-section 
for the photon-photon production of Be- 

At the LHC energies, the parton subprocess of gluon-gluon fusion gg — >■ 

+ b + c dominates in the hadron-hadron production of Be mesons. In the 
leading approximation of QCD perturbation theory it requires the calculation 
of 36 diagrams in the fourth order over the coupling constant. In this case, 
there are no isolated gauge-invariant groups of diagrams, which would allow 
the interpretation similar to the consideration of Be production in e^e~- 
annihilation and photon-photon collisions. 

The calculation for the complete set of diagrams of the 0(af)-contribution 
allows one to determine a value of the transverse momentum p^*", being the 
low boundary of the region, where the subprocess of gluon-gluon i?c-nieson 
production enters the regime of factorization for the hard production 66-pair 
and the subsequent fragmentation of 5-quark into the bound (6c)-state, as it 
follows from the theorem on the factorization of the hard processes in the 
perturbative QCD ||T0[| . 

The p™" value turns out to be much greater than the Mb^ mass, so 
that the dominant contribution into the total cross-section of gluon-gluon 
-Bc-pro duct ion is given by the diagrams of nonfragmentational type, i.e. by 
the recombination of heavy quarks. Furthermore, the convolution of the 
parton cross-section with the gluon distributions inside the initial hadrons 
leads to the suppression of contributions at large transverse momenta as 
well as the subprocesses with large energy in the system of parton mass 
centre, so that the main contribution into the total cross-section of hadronic 
i?c-production is given by the region of energies less or comparable to the 
-Bc-meson mass, where the fragmentation model can not be applied by its 
construction. Therefore, one must perform the calculations with the account 
for all contributions in the given order under consideration in the region near 
the threshold. 

In this work we calculate the both total and differential cross-sections for 
the production of 5*- and P-wave states of Be mesons in the framework of 
leading order of the QCD perturbation theory. Below we give a quite com- 
plete set of the Be characteristics. The basic spectroscopic characteristics of 
the (be) family and the branching ratios of different decay modes of the pseu- 
doscalar ground state, which can be used for the extraction of the Be signal 
from the hadronic background, are presented in section 1. The calculations 
of characteristics in the hadronic Be production are analyzed in section 2. In 
Conclusion the obtained results are summarized. 
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1 The spectroscopy and decays of Be 



1.1 Level masses and coupling constants 

As a system containing two heavy quarks, the (6c)-mesons can be rehably de- 
scribed in the framework of phenomenological nonrelativistic potential mod- 
els [|rT] and in QCD sum rules |T^, which are based on (i) the formalism of 
the expansion over the inverse heavy quarks mass, due to a small ratio of the 
typical quark virtualities, being of the order of the confinement scale, to the 
mass, and (ii) the expansion over a small relative velocity of the quark motion 
in the bound state. The perturbative QCD allows one to take into account 
hard gluon corrections to the currents of effective quarks, defining the lead- 
ing approximation in these expansions ]13|]. This way of consideration leads 
to a quite precise and reliable description of the families of (cc)-charmonium 
and (66)-bottomonium, whose spectroscopic properties are in detail studied 
experimentally P |. 



The effective gluon potential of static quarks does not depend on the 
flavour of heavy quark. The different phenomenological models lead to the 
potential, whose form in the range of average distances between the quarks 
for ground and excited levels: 0.1 < r < 1 fm, really does not depend on 



the model |jT5|]. The experimental data on the (cc) and (bb) systems point 
out the energy level regularity expressed in the approximate quark-flavour 
independence of the difference between energies of excitations. According to 
the Hell-Mann - Feynman theorem, this regularity is caused by the constant 
value of average kinetic energy of the heavy quarks. This energy does not 
depend on the quark flavours and number of excitation. The independence 
becomes exact for the logarithmic potential |jl6| . The weak dependence of the 



difference between energy levels can be taken into account while considering 
Martin potential: V{r) = v4(r/ro)" + C 0]- 

The mass values for the different states of (6c) system in Martin potential 
with account for the spin-dependent splitting of riL-levels are presented in 
table 1. The splittings are calculated in the second order of l/mg-expansion 
over the inverse heavy quark mass with the account for the vector coupling 
of the effective one-gluon exchange between the quarks and scalar coupling 



for the part of potential confining the quarks in the bound state [0. The 
effective value of the gluon coupling constant at the scale of the average 
momentum transfer between the quarks can be determined over the known 
experimental value for the splitting of charmonium IS'-level with the taking 
into account the renormalization group evolution towards the {be) energies 

i- 

Thus, in the (be) system below the threshold of decay into the heavy 
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Table 1: The masses of bound (6c)-states below the threshold of decay into 
the pair of heavy mesons BD, in GeV, in the model with Martin potential 
and the ST-potential motivated by QCD with the account for the two-loop 
evolution of the coupling constant at short distances [19]. The spectroscopic 
notations of states are n^^^Lj, where jc is the total angular momentum of 
c-quark, n is the principal quantum number, L is the orbital angular mo- 
mentum, J is the total spin of meson. 



state 


Martin 


i 




BT 








6.253 


6.264 




6.317 


6.337 


2'So 


6.867 


6.856 


2^ Si 


6.902 


6.899 


2'Po 


6.683 


6.700 


2P 1+ 


6.717 


6.730 


2P 1'+ 


6.729 


6.736 


2^P2 


6.743 


6.747 




7.088 


7.108 


3P 1+ 


7.113 


7.135 


3P 1'+ 


7.124 


7.142 


3'P2 


7.134 


7.153 


3D 2- 


7.001 


7.009 




7.007 


7.005 


3^Di 


7.008 


7.012 


3D 2'- 


7.016 


7.012 



meson BD pair, one can expect the presence of 16 narrow states of the 
IS*, 2S,2P, 3P and 3P'-levels, which cascadely transform into the ground 
pseudoscalar state with the mass m^^ = 6.25 ±0.03 GeV due to the emission 
of photons and vr-meson pairs, since the annihilation (6c) decays can occur 
only due to the weak interaction and, hence, are suppressed for the excited 
levels. The total widths of the excited states and the branching ratios for 
the radiation modes of decay are shown in table 2. 

The potentials motivated by QCD have the linear growth at large dis- 
tances and the coulumb-like behaviour at short ones. The form of the po- 
tential in the region, where the perturbative regime changes into the nonper- 
turbative one with the increase of the distance between the quarks, coincides 
with the form of the logarithmic or power potentials, so that the accuracy 
for the prediction of the energy levels in the heavy quarkonia, particular in 
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the (be) system, is determined by the value of 30 MeV. 

The global properties of these potentials, i.e. their asymptotic behaviour 
in the bound points of r — oo, r — 0, are significant for the determination of 
the coupling constants of the states, such as the leptonic coupling constants 
/ for the nS'-levels, for example. In the leading approximation, the / value 
does not depend on the spin state of the level and it is determined by the 
value of the radial wave function at the origin, -R(O), 



Taking into account the hard gluon corrections [^, |2T|, the constants of 
vector and pseudoscalar states equal 



/y-Aii+^r^^in^-*"'")), (1) 

V vr Vmi+m2 m2 J J 

where m\^2 are the quark masses, 5^ = 8/3, = 2, and the QCD coupling 
constant is determined at the scale of the gluon virtuality, given by the quark 
masses function . 



The values of the heavy quarkonium wave functions are significant in the 
consideration of production and decays of their states in the framework of 
quark models of the mesons. They are shown in table 3 for S- and P-wave 
states in the model with Martin potential. However, as it was mentioned 
above, the accuracy of potential models for the values under consideration is 
quite low, and the uncertainty is expressed by a factor of two. 

The QCD sum rules allow one to determine the leptonic constants for the 
heavy quarkonium states with a much better accuracy. Standard schemes of 
the sum rules give an opportunity to calculate the ground state constants 
for vector and pseudoscalar currents with the account for corrections from 
the quark-gluon condensates, which have the power form over the inverse 



heavy quark mass [|T2| . There is a region of the momentum numbers for the 
spectral density of the two-point current correlator, where the condensate 
contributions are not significant. In this region, the integral representation 
for the contribution of resonances lying bellow the threshold of decay into the 
pair of heavy mesons, allows one to use the regularity of the quarkonium state 
density mentioned above, and to derive the scaling relations for the leptonic 
constants of the ground state quarkonia with different quark contents and 



for the excited states 211. Thus, for vector states we have 



Mn \mJ V 4mi2 y n ' ^ ^ 
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Table 2: The total widths of excited bound (6c)-states below the threshold 
of decay into the meson BD-pair in the model with Martin potential and the 
branching ratios of the dominant decay modes. 



state 


Ttot, KeV 


dominant decay 
mode 


BR, % 


l^Si 


0.06 


l'5o + 7 


100 




67.8 


l^^o + rni 


74 


2^Si 


86.3 


l^Si + nir 


58 




65.3 


1^5i+7 


100 


2P 1+ 


89.4 


1^5i+7 


87 


2P 1'+ 


139.2 


1^^0 + 7 


94 




102.9 


1^^1+7 


100 




44.8 


2^^1+7 


57 


3P 1+ 


65.3 


2^51+7 


49 


3P 1'+ 


92.8 


2^^0 + 7 


63 




71.6 


2^5i+7 


69 


3D 2- 


95.0 


2P 1+ + 7 


47 




107.9 


23P2 + 7 


71 




155.4 


2^Po + 7 


51 


3D 2'- 


122.0 


2P 1'+ + 7 


38 



where mi2 = 1711^12/ {rrii + 1112) is the reduced mass of quarks, and the con- 
stant c is determined by the average kinetic energy of quarks, the QCD 
coupling constant at the scale of average momentum transfers in the system 
and the hard gluon correction factor to the vector current. Numerically, in 
the method accuracy the c value turns out to depend on no quark flavour 
and excitation number in the system. Relation (Q) is in a good agreement 
with the data on the coupling constants for the families of ip and T particles 
and, thus, it can be a reliable basis for the prediction of leptonic constants 
for the Pc-nieson family, as it is shown in table 4, where one can also find 
the constants for the pseudoscalar states of nS'-levels, which are determined 
by the relation (see pT[] ) 

2as\ mi + 1112 
~3ti) Wn ■ 

The P(0) values, being the parameters of static quark models in the 
consideration of decays and production of the mesons, can be determined over 
the leptonic constants calculated in the sum rules. In this way, the additional 
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Table 3: The wave function characteristics, Rns{0) (in GeV^/^) and -R^p(O) 
(in GeV^/^), obtained from Schrodinger equation with Martin potential and 
BT. 



n 


Martin 


i 


BT 


i 


Risio) 

R2S{0) 
^2P(0) 
^3P(0) 


1.31 
0.97 
0.55 
0.57 


1.28 
0.99 
0.45 
0.51 



Table 4: The leptonic constants for the vector and pseudoscalar states of 
raS'-levels in the (6c)-system, /„ and f^, calculated in the sum rules resulting 
in the scaling relation, and Rns{0), calculated over /„. The accuracy is equal 
to 6%. 



n 


fn, MeV 


f^, MeV 


i?„5(0), GeV'r2 


1 


385 


397 


1.20 


2 


260 


245 


0.85 



uncertainty appears. It is related with a choice of the scale fixing the QCD 
coupling constant in the factor accounting for the hard gluon correction. The 
dependence on this choice points out the significance of high orders in the 
perturbative approximation. Assuming that in the MS-scheme the scale is 
close to /i^ = mim2e~^^/^^ |22[, one obtains the values of -R(O) shown in table 



4, where one can see that the predictions of sum rules and the potential 
models of Martin and BT are in a reasonable agreement with each other. 

As for the states lying above the threshold of decay into the heavy meson 
BD pair, the width of B*'^{3S) B~^D^, for example, can be calculated 
in the framework of sum rules for the meson currents, where the scaling 
relation takes place for the g constants of similar decays of T(4S') B^B~ 
and ^/^(3770) ^ D+D- [El, 



const. 



r_ /4mi2\ 

m\m) 



The relation is caused by the dependence of energy gap between the vector 
and pseudoscalar heavy meson states: AMi^2 ■ Mi^2 = const., where Mi^2 are 
the meson masses in the final state, mi2 is their reduced mass. The width of 
this decay has a strong dependence on the B*^{3S) mass, and at M = 7.25 
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GeV it is equal to F = 90 ± 35 MeV, where the uncertainty is determined by 
the accuracy of method used. 

1.2 The Be meson decays 
1.2.1 Life-time 

The 5c-nieson decay processes can be subdivided into three classes: 1) the 
6-quark decay with the spectator c-quark, 2) the c-quark decay with the 
spectator 6-quark and 3) the annihilation channel — >• l~^i'i{cs,us), where 
I = e, /i, T. In the b ccs decays one separates also the Pauli interfer- 
ence with the c-quark from the initial state. In accordance with the given 
classification, the total width is the sum over the partial widths 

r(5e ^ X) = T{b ^X)+ r(c ^ X) + r(ann.) + r(PI) . (3) 

For the annihilation channel the r(ann.) width can be reliably estimated in 
the framework of inclusive approach, where one take the sum of the leptonic 
and quark decay modes with account for the hard gluon corrections to the 
effective four-quark interaction of weak currents. These corrections result in 
the factor of ai = 1.22±0.04 The width is expressed through the leptonic 
constant of Jb^ = fi ^ 400 MeV. This estimate of the quark-contribution 
does not depend on a hadronization model, since the large energy release of 
the order of the meson mass takes place. From the following expression one 
can see that one can neglect the contribution by light leptons and quarks, 

F(ann.) = f^lHd VI Mm,2(l - m^/mlf ■ Q , (4) 

i=T,C 

where Cr = 1 for the r"''z/T--channel and C = 3\Vcs\'^a1 for the cs-channel. 

As for the nonannihilation decays, in the approach of the operator product 
expansion for the quark currents of weak decays one takes into account 
the a^-corrections to the free quark decays and uses the quark-hadron duality 
for the final states. Then one considers the matrix element for the transition 
operator over the bound meson state. The latter allows one also to take 
into account effects caused by the motion and virtuality of decaying quark 
inside the meson because of the interaction with the spectator. In this way 
the b — >• CCS decay mode turns to be suppressed almost completely due to 
the Pauli interference with the charm quark from the initial state. Besides, 
the c-quark decays with the spectator 6-quark is essentially suppressed in 
comparison with the free quark decays because of a large bound energy in 
the initial state. 
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Table 5: The branching ratios of the Be decay modes calculated in the frame- 
work of inclusive approach and in the exclusive quark model with the pa- 
rameters \Vbc\ = 0.040, fe, = 0.47 GeV, = 0.54 GeVJe, = 0.3 GeV, 
nib = 4.8 — 4.9 GeV, rUc = 1.5 — 1.6 GeV, = 0.55 GeV. The accuracy is 
about 10%. 



Be decay mode 


Indus., % 


Exclus., % 


b cl^i'i 


3.9 


3.7 


b — > cud 


16.2 


16.7 




25.0 


25.0 


c — > sl^ui 


8.5 


10.1 


c —>■ sud 


47.3 


45.4 


Ec^s 


64.3 


65.6 


B+ r+Ur 


2.9 


2.0 


Bt ^ cs 


7.2 


7.2 



In the framework of exclusive approach it is necessary to sum widths 
of different decay modes calculated in the potential models 0, While 
considering the semileptonic decays due to the b — > cl^vi and c sl^Pi 
transitions one finds that in the former decays the hadronic final state is 
practically saturated by the lightest bound IS'-state in the (cc)-system, i.e. 
by the r]c and J/ip particles, and in the latter decays the IS'-states in the 
(6s)-system, i.e. Bs and -B*, only can enter the accessible energetic gap. The 
energy release in the latter transition is low in comparison with the meson 
masses, and, therefore, a visible deviation from the picture of quark-hadron 
duality is possible. Numerical estimates show that the value of B^ — > {bs)l'^i>i 
decay width is two times less in the exclusive approach than in the inclusive 
method, though this fact can be caused by the choice of narrow wave package 
for the 5^*-* mesons in the quark model, so that /^^ ~ 150 MeV, while in the 
limit of static heavy quark, one should expect a larger value for the leptonic 
constant |2^. This increase will lead to the widening of the wave package 
and, hence, to the increase of the overlapping integral for the wave functions 
of and 5^*) (see table 5). 

Further, the b — >■ cud channel, for example, can be calculated through the 
given decay width of 6 — cl'^vi with account for the color factor and hard 
gluon corrections to the four-quark interaction. It can be also obtained as a 
sum over the widths of decays with the (M(i)-system bound states. 

The results of calculation for the total B^. width in the inclusive and 
exclusive approaches give the values consistent with each other, if one takes 
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into account the most significant uncertainty related with the choice of quark 
masses (especially for the charm quark), so that finally we have 

r(5+) = 0.55 ±0.15 psec. (5) 
1.2.2 Exclusive decays 

The consideration of exclusive i?c-decay modes supposes an introduction of 
model for the hadronization of quarks into the mesons with the given quan- 
tum numbers. The QCD sum rules for the three-point correlators of quark 
currents ^ and potential models 0, ^] are among of those hadronization 



models. A feature of the sum rule application to the mesons containing two 
heavy quarks, is the account for the significant role of the coulumb-like ag/v- 
corrections due to the gluon exchange between the quarks composing the me- 
son and moving with the relative velocity v. So, in the semileptonic decays of 

— > ilj{rjc)l^Vh the heavy {Q1Q2) quarkonium is present in the both initial 
and final states, and, therefore, the contribution of coulumb-like corrections 
exhibits in a specially strong form. The use of tree approximation for the 
perturbative contribution into the three-point correlator of quark currents 

leads to a large deviation between the values of transition form-factors, 
calculated in the sum rules and potential models, respectively. The account 
for the a^/'i'-corrections removes this contradiction Thus, the meson 



potential models, based on the covariant expression for the form-factors of 
weak Be decays through the overlapping of quarkonium wave functions in the 
initial and final states, and the QCD sum rules give the consistent description 
of semileptonic Sc-meson decays. 

Further, the hadronic decay widths can be obtained on the basis of as- 
sumption on the factorization of the weak transition between the quarkonia 
and the hadronization of products of the virtual *"^-boson decay . The 
accuracy of factorization has to raise with the increase of VT-boson virtuality. 
This fact is caused by the suppression of interaction in the final state. In this 
way, the hadronic decays can be calculated due to the use of form-factors for 
the semileptonic transitions with the relevant description of W* transition 
into the hadronic state. 

Let us consider the amplitude of 5+ — * Mxe'^i'e transition with the weak 
decay of quark 1 into the quark 2 

A=^V,,l,H^, (6) 

where Gi? is the Fermi constant, V12 is the Kobayashi-Maskawa matrix ele- 
ment. The leptonic current is determined by the expression 

l^l = e(gi)7/.(l - 75)^^(92) , (7) 
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where qi^2 are the momenta of lepton and neutrino, correspondingly; (gi + 
^2)^ = t. The value in is the matrix element for the hadronic 
current 

J^ = V;-A^ = Qi7^(l- 75)^2. (8) 

The matrix element for the decay into the pseudoscalar P state can be 
written down in the form 

< B,{p)\A,\Pik) >= + k), + - k), , (9) 

and for the transition into the vector V meson with the My mass and A 
polarization one can express the matrix element in the form 



<B,{p)\J^\Vik,\)> = -{M + Mv)Aiit) e 



(A) 



Relations ( [T0| ) define the form-factors for the B^ Mxe^Ve transitions, 
so that for the massless leptons, F_ and ^3 do not contribute into the matrix 
element in (^. In the covariant quarkonium model with the wave functions 
of oscillator type @] one can easily obtain 



F,(t) = l(,„,H-™,)^^i-{p(t). (11) 

F-{t) = -i(mi-m2 + 2m.p)^/^i-?p((), (12) 

where rusp is the mass of spectator quark, and the ^(t) function has the form 

where Mx is the meson mass, and lOx is the wave function parameter of the 
recoil meson, 
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Table 6: The branching ratios of exclusive Be decay modes, calculated in 
the framework of covariant quark model with the parameters |Vfec| = 0.040, 
/b, = 0.47 GeV, = 0.54 GeV, fe, = 0.3 GeV, rub = 4.8 - 4.9 GeV, 
rric = 1.5 — 1.6 GeV, = 0.55 GeV. The accuracy equals 10%. 



Be decay mode 


BR, % 


Be decay mode 


BR, % 






2.5 




1.2 


B*J+ui 




6.2 


Bsl+ui 


3.9 






0.2 




0.2 






5.2 


Bsn+ 


5.5 


'ipp+ 




0.6 


VcP^ 


0.5 






22.9 


BsP+ 


11.8 



is the maximum square of the leptonic pair mass. For the vector state 
Mv = Mx we get 



V[t) ^ \iM + My)^^Uv{t). (15) 

Mt) = \{M + My){l-2mjM)sl^^J,y{t), (17) 

A,{t) = -UM + My){l + 2mjM)\^—Cv{t). (18) 

The branching ratios of some exclusive i?c-decay modes are shown in table 
6. 

A decrease of the invariant mass for the hadron system results in an 
increase of the recoil meson momentum. This causes the problem of appli- 
cability for the formalism of overlapping for the quarkonium wave functions, 
because, in this kinematics, the narrow wave packages are displaced relative 
to each other in the momentum space into the range of distribution tails. In 
this situation one has to take into account a hard gluon exchange between 
the quarkonium constituents, which destroys the spectator picture of weak 
transition in the potential approach. For example, the widths of B^ —>■ ipir^ 
and B'^ ^c^^^ decays have the following forms |2^ 

r(fi, -.i^n) - G,\Vbe\ ^^fJB.U[jf^) (M-m^)2^2(i?>, 
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Table 7: The branching ratios of exclusive Be decay modes at the large 
momenta of recoil meson in the framework of hard gluon exchange approxi- 
mation. The uncertainty equals 50%. 



Be decay mode 


BR., % 


Be decay mode 


BR, % 




2.0 




0.7 




5.6 


VcP^ 


2.0 



-u (5-m2/M2)2 



As one can see in equations (p!9D-(P0|), the basic uncertainty in the width 
estimation is connected with the scale choice for the "running" QCD coupling 
constant, which supposes the renormalization of the six-quark operator for 
the transition between the quarkonia in the next-to-tree approximation. This 
problem requires a carrying out a large volume of cumbersome calculations. 
However, even in this stage of evaluations one can stat that the corrections 
with "large" logarithms expect to appear from the contributions related with 
the renormalization of four-quark weak interaction (the Oi factor) and in the 
gluon propagator standing in the exchange between the quarks. The latter 
leads to the fixing of scale at the gluon virtuality. Numerically, this 
virtuality is determined by the recoil meson momentum, and it is slightly 
greater than 1 GeV, so that one can talk on a quite hard gluon exchange, 
which gives the values of branching ratios of approximately one order of 
magnitude greater than in the potential approach, as one can find in table 7. 
However, the conservative estimate of uncertainty in the hard approximation 
is close to 50%. 

As for the extraction of Be signal in the hadronic background, the decay 
modes with ip in the final state are the most preferable, because the latter 
particle can be easily identified by its leptonic decay mode. This advantage 
is absent in the Be decay modes with the final state containing the r]e or B^*^ 
mesons, which are the objects, whose experimental registration is impeded by 
a large hadron background. From the values of branching ratios shown above 
one can easily obtain, that the total probability of ip yield in the Be decays 
equals BR(i?^ ^-^) = 0.24. It is worth to note that the key role in the 
Be signal observation plays the presence of the vertex detector, which allows 
one to extract events with the weak decays of long-lived particles containing 
heavy quarks. In the case under consideration it gives the possibility to 
suppress the background from the direct ip production. In the semileptonic 
B^ — > ipl^ui decays the presence of vertex detector and large statistics of 
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events allows one to determine the i?c-nieson mass and to separate the events 
with its decays from the ordinary i?+-meson decays, which have no ipl^ mode. 
In the ipn^ decay the direct measurement of Be mass is possible. The detector 
efficiency in the reconstruction of three-particle secondary vertex (/+/~ from 
decays of ijj and 7r+ or /"*") becomes the most important characteristics here. 
The low efficiency of the LEP detectors (e ~ 0.15), for example, makes the Be 
observation to be probably unreachable in the experiments at the electron- 
positron collider pO . 



2 The characteristics of production at LHC 

As was mentioned above, the consideration of mechanisms for the hadronic 
production of different spin i?c-states is based on the factorization of hard 
parton production of heavy quarks (bbcc) and soft coupling of (be) bound 
state [|^. In the first stage of description, the hard subprocess can be reli- 
ably calculated in the framework of QCD perturbation theory, while in the 
second stage the quark binding in the heavy quarkonium can be described in 
the nonrelativistic potential model assigned to the (6c)-pair rest system. The 
latter means that one performs the integration of the final quark state over 
the quarkonium wave function in the momentum space. Since the relative 
quark velocity inside the meson is close to zero, the perturbative matrix ele- 
ment can be expanded in series over the relative quark momentum, which is 
low in comparison with the quark masses determining the scale of virtualities 
and energies in the matrix element. In the leading approximation one con- 
siders only the first nonzero term of such expansion, so that for the S-wave 
states the matrix element of the parton subprocess for the Be production is 
expressed through the perturbative matrix element for the production of four 
heavy quarks {gg — >■ bbcc) with the corresponding projection to the vector or 
pseudoscalar spin state of (6c)-system, being the color singlet, and through 
the factor of radial wave function at the origin, Rns{0), for the given quarko- 
nium. The perturbative matrix element is calculated for the [be) state, where 
the quarks move with the same velocity, i.e. one neglects the relative motion 
of b and c. 

For the P-wave states, the potential model gives the factor in the form 
of first derivative of the radial wave function at the origin, i?^^(0). In the 
perturbative part, one has to calculate the first derivative of the matrix 
element over the relative quark momentum at the point, where the velocities 
of quarks, entering the quarkonium, equal each to other. 

Thus, in addition to the heavy quark masses, the values of ^^^(O), -R^2.(0) 
and as are the parameters of calculation for the partonic production of Be 
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meson. In calculations we use the wave function parameters equal to the 
values shown in table 4 and i?2p(0) = 0.50 GeV^^^. The value of -R(O), as 
was shown above, can be related with the leptonic constant, /, so that we 
have 

hs = 0.47 GeV, hs = 0.32 GeV. 

At large transverse momenta of the Be meson, pt ^ Mb^, the production 
mechanism enters the regime of 6-quark fragmentation (see fig. 1), so that 
the scale determining the QCD coupling constant in hard bb production, 
is given by nf^^ ~ M^^ + p^, and in the hard fragmentation production of 
the additional pair of heavy quarks cc we get /xgc ~ This scale choice is 
caused by the high order corrections of perturbation theory to the hard gluon 
propagators, where the summing of logarithms over the virtualities leads to 
the pointed /i values. Therefore, the normalization of matrix element is 
determined by the value of as{fibb)'^s{fJ'cc) ~ 0.18 • 0.28. In calculations we 
use the single combined value of as = 0.22. 

In the approximation of the weak quark binding inside the meson one 
has Mb^ = rrib + rric, so that the performable phase space in calculations is 
close to physical one at the choices of rrib = 4.8 GeV, rric = 1.5 GeV for 
the IS'-state, rrib = 5.1 GeV, rric = 1-8 GeV for the 2S'-state, rrib = 5.0 
GeV, rric = 1-7 GeV for the 2P-state. The value of the hadron-hadron cross- 
section is determined by the convolution of parton distributions with the 
cross-sections of parton subprocesses. At the LHC energies the subprocess of 
gg — > B^bc is dominant in the production of Be meson because of the more 
high gluon-gluon luminosity. The total cross-sections of the subprocesses 
for the different spin states of Be meson are calculated numerically, and 
their dependence on the subprocess energy, V^, can be written down in the 
following analytic approximation 

a{gg ^ Bc{lS)bc) = 550.0 1 - — I K cj\ 



a(,, ft(2P) ^ 25.0 (l - ( " p^?2, 



in the region of 2{mb + rric) < < 400 GeV. 

The using of CTEQ4M parameterization for the structure functions of 
nucleon |^ leads to the total hadronic cross- sections for the Be mesons. 



as shown in table 8. After the summing over the different spin states, the 
total cross-sections for the production of P-wave levels is equal to 7% of the 
S'-state cross-section. 

In fig. 2 the da/dpT distributions over the transverse momentum for the 
yields of IS-states of Be and B* are shown at the 8 and 16 TeV energies. The 
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Table 8: The total cross-sections for the hadronic production of different spin 
states of the Be mesons at the LHC energies. 



nLj 


a (8 TeV), nb 


a (16 TeV), nb 


ISo 


21.961 


46.104 


ISi 


55.039 


115.493 


2So 


4.818 


10.115 


2Si 


12.022 


25.226 



maximum in the distribution is reached at pt ~ Mb^- Prom the distribution 
over rapidity (fig. 3) one can sec that the production of Be meson takes place 
basically in the central region. The topology of events with the Bc-raeson 
production somewhat differs from the configuration for the B meson. In the 
system of mass centre for the final Be + b + c state, the Be and b particles 
mainly move in opposite directions. Otherwise, as it is shown in fig. 4, the 
c-quark moves in the same direction as the Be meson. The distribution over 
the angle between the Be meson and c-quark has the sharp maximum in the 
region of 0.95 < cos^^ < 1. 

Thus, the production of Be meson is associated by the production of D 
meson moving in the same direction as Be, i.e. the D velocity is almost 
parallel to the Be one, 9 < 18°. The average transverse momentum of D 
meson is less than that of Be (see fig. 5) in this event . As the registration 
of Be meson supposes the observation of secondary vertex, two vertices must 
be searched for in the same jet in the events with Be- These are the vertices 
from the Be and D meson decays. The correlation between the energy of 
Be meson and c-quark associated with him is shown in fig. 6. A rough 
estimate shows that the c-quark energy is approximately the same as that of 
Be- Taking into account the fact, that after the fragmentation the average 
D-meson energy equals (Ed) ~ 0.6{Ee), and also the average lifetime of Be 
is approximately two times less than that of the D'^ meson it is easily to 
obtain the estimates of the distances between the primary interaction point 
and the decay vertices of D meson and Be- 

Ec rx c 

xd^ YF~ 0-6crD, xb, ^ -ji— ctb, , 

Mb^ 

i.e. 

~ 3.6 . 

Thus, in the events with Be one should expect also, for example, the asso- 



17 



ciated production of mesons, having approximately three times greater 
displacement from the primary vertex of interaction for the particles from the 
beams, than the displacement of Be, and the direction of D-meson motion 
approximately coincides with that of Be- This configuration can be useful in 
the extraction of signal from the hadronic background. 

Conclusion 

In this paper we have considered the production mechanism for the different 
spin S- and P-wave states of the Be mesons in the hadron interactions at the 
LHC energies on the basis of calculation for the complete set of diagrams for 
the leading 0(Q;^)-contribution of the perturbative QCD for the hard pro- 
duction of heavy quarks and in the potential model of soft binding of the 
b and c quarks into the meson. One has found that the regime of 6-quark 
fragmentation into the Be meson is delayed to the region of large transverse 
momenta, pr > QMb^ ^ Mb^, so that the diagrams of the nonfragmen- 
tational type, i.e. the recombinations, give the dominant contribution into 
the total cross-section for the hadronic production of Be- At the maximum 
luminosity £ = 10^"^ cm~^c~^, being planned, one should expect the yield of 
Be at the level of 3 • lO-^*^ events per year at ^/s — 8 TeV and 6 • 10^° events 
per year at a/s — 16 TeV. 

The theoretical predictions of spectroscopic characteristics for the Bc- 
meson family and their decays allow one to carry out an object-oriented 
search for the events with Be and to extract them from the hadron back- 
ground. For example, in the semileptonic decay mode of B^ — > ipl'^i'i with 
the account for the branching ratio for the leptonic decay of ■0 — I'^l' and 
the efficiency of the registration for the secondary vertex with three leptons 
e = 0.1, one should expect 12 ■ 10^ reconstructed events with the Be mesons 
at y/s = 16 TeV. Even in the regime of low energies {^/s = 8 TeV) and lu- 
minosity (£ = 10^^ cm~^c~^) it is possible to reconstruct 6 • 10^ semileptonic 
Be decays only in the muon channel, for example. The probability of the 
B^ —>■ %l)TT^ decay has the uncertainty related with the large momentum of 
the recoil ijj meson. In the regime of low energies and luminosity the number 
of the reliably detected events is about (0.4 — 4) • 10*^ at the same efficiency 
of registration of the secondary vertex (e = 0.1). It is worth to mention that 
we have pointed out the number of events with no account for the detector 
cuts off the angles and momenta in conditions of particular facilities at LHC. 
Nevertheless, the search for Be at LHC seems to be the most promising for a 
positive result, since, for example, at the LEP facilities the observation of Be 
is extremely problematic due to the low yield of Be- At HERA-B a strong 
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threshold effect for the hadron production of Be will appear because one has 
to produce the additional pair of heavy B and D mesons in the final state. 
So, the statistics at HERA-B will be at the same level as that of LEP, and, in 
addition, it will have conditions of a more intensive hadron background. A 
low ratio of the signal to a background also does not yet allow the CDF col- 
laboration (FNAL) to observe the events with Be because of a low efficiency 



of the reconstruction for the secondary vertex [32 



Thus, the search for the Be meson at LHC thinks to be a quite interesting 
and solvable problem. 

This work is in part supported by Russian Foundation for Basic Research, 
grant 96-02-18216. 
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Figure 1: The differential cross-section for tlie i?^*-* meson production in 
gluon-gluon collisions as calculated in the perturbative QCD over the com- 
plete set of diagrams in the 0{af) order at 200 GeV. The dashed and solid 

histograms present the pseudoscalar and vector states, respectively, in com- 
parison with the results of fragmentation model shown by the corresponding 
smooth curves. 
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Figure 2: The distributions over the transverse momentum of B* (sohd hne) 
and Be (dashed hne) in proton-proton interactions at the energies of 8 TeV 
(a) and 16 TeV (b). 
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Figure 3: The distributions over the rapidities of B* (sohd hne) and Be 
(dashed hne) in proton-proton interactions at the energies of 8 TeV (a) and 
16 TeV (b). 
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Figure 4: The da/dcos6 distribution over the angle of 6 between the 
meson and c-quark in the system of mass centre for the coUiding proton 
beams at the energy of 8 TeV. 
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Figure 5: The distribution over the transverse momentum of the c-quark 
associated with the B,, meson in proton-proton interactions at the energy of 
8 TeV. 
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Figure 6: The correlation between the energies of the Be meson and c-quark 
associated with the former at the LHC energy of 8 TeV. 
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